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Abstract. In this paper we report multi-frequency measurements of pulse broadening times (ra) for nine medium 
dispersion measure (DM ~ 150 — 400 pc cm -3 ) pulsars observed over a wide frequency range. The low frequency 
data at 243, 325 and 610 MHz are new observations done with the Giant Metrewave Radio Telescope (GMRT). 
The frequency dependence of t<j for all but one (PSR B1933+16) of our sources is consistent with the Kolmogorov 
spectrum of electron density fluctuations in a turbulent medium. PSR B1933+16, however, shows a very flat 
spectrum as previously observed for high DM pulsars. Our observations combined with earlier published results 
enable us to study the spectral index of t^ over the whole observed DM range. While the spectral properties are 
generally consistent with a Kolmogorov spectrum, pulsars seen along line-of-sights towards the inner Galaxy or 
complex regions often show deviations from this expected behaviour. 
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1. Introduction 

The free electron density distribution in the interstellar 
medium (ISM) can be decomposed into three different 
regimes: the spiral arms consisting of gaint HII regions, 
the inner disk comprising of a dense ionized region and the 
thick disk filled with diffuse electron gas (Taylor & Cordes 
1993). The pulsar signal traversing the ISM suffers inter- 
stellar dispersion quantified in terms of dispersion mea- 
sure, DM, which is the integrated electron column density 
towards the pulsar (at distance D), i.e. DM = f Q n dl. 
Fluctuations of the electron density along the line of sight 
(LOS) give rise to several observable scattering effects, 
which manifest themselves as observational properties like 
angular broadening, temporal pulse broadening and scin- 
tillation of pulsars (see Rickett 1990 for a review). 

Scattering causes propagation of signals along a variety 
of different ray paths with different geometrical lengths, 
so that a pulse, which has left the pulsar at one instant, 
arrives at the observer over a finite time interval. For a 
gaussian distribution of irregularities and applying the 
thin screen approximation, the pulse broadening function 
(PBF) of the ISM can be well described by an exponential 



decay of the pulse, i.e. PBF(t) = cxp(— t/r d ), where Td is 
called the pulse broadening time. (Scheuer 1968). Random 
interference among the different paths produces a diffrac- 
tion pattern in the plane of the observer. This pattern 
decorrelates over a characteristic bandwidth Ais^- Both 
Td and Av<i strongly depend on frequency v (i.e. Td oc v~ a 
and A^d oc v a ) and are related to each other as 



27T Td Al/d = C\ 



(1) 



The constant C\ is of the order unity but changes for 
different models and geometries for the electron density 
fluctuations (e.g. Lambert & Rickett 1999). However, until 
now there exists no direct measurement of Td and A^d at 
the same observing frequency for obtaining C± . 

The strength of scattering of radio waves caused by 
electron density fluctuations in the ISM has been subject 
of detailed studies since the discovery of scintillation in 
pulsar signals (Scheuer 1968). A commonly used descrip- 
tion of the scattering strength is to attribute a power law 
to the electron density spectrum with a large range be- 
tween "inner" and "outer" scales, k, 

r<2 / 

( " II 



and k„ ; i.e. 



p*M 



exp 



(2) 
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where q is the magnitude of the three-dimensional wave- 
number and C^ is the fluctuation strength for a given 
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LOS (e.g. Rickett 1977). For k < q < k t , one ob- 
tains a simple power-law model with a spectral index /?, 
i.e., P ne (q) = C% e q-P and a = 2/3/(0 - 2). For a pure 
Kolmogorov spectrum of density irregularities, f3 = 11/3, 
we expect a — 4.4. Attempts to reconstruct the electron 
density spectrum from various observations (Armstrong 
et al. 1995) have led to the suggestion that over 5 decades 
of wavenumber range 10 -13 m" 1 < q < 10~ 8 m" 1 the 
spectrum can be approximated by a power law with an 
average index ( (3 ) ~ 3.7, close to the value expected from 
the Kolmogorov process. However, there is some evidence 
for deviations from this picture, over some parts of the 
wavenumber range (see, for example, Gupta 2000 and ref- 
erence therein). 

Furthermore, it is uncertain how the scattering ma- 
terial and processes operate in different directions of 
the Galaxy. Multi-frequency scatter broadening measure- 
ments have been used to determine a for several LOSs in 
the Galaxy. For LOSs with low DM's (< 50 pc cm -3 ) there 
is typically good agreement with the Kolmogorov spec- 
trum (Cordes et al. 1985; Johnston et al. 1998). This is also 
supported by estimates of a in the local ISM using diffrac- 
tive and refractive interstellar scintillation (ISS) observa- 
tions of pulsars (Bhat et al. 1999). These local LOSs are 
most probably sampling the diffuse, homogeneous electron 
density component of the Galaxy where the Kolmogorov 
process hold true. High DM (> 400 pc cm -3 ) pulsars to- 
wards the inner parts of the Galaxy, in contrast, have an 
average (a) = 3.44 ±0.13, showing significant departures 
from the Kolmogorov spectrum (Lohmer et al. 2001, here- 
after Paper I). We interpreted this phenomenon in Paper 
I as anomalous scattering as discussed by Cordes & Lazio 
(2001). In this picture the observed anomalous behaviour 
can be explained by invoking scattering caused by multiple 
scattering screens with anisotropic irregularities and finite 
transverse extent along the LOS. As a consequence, less 
radiation reaches the observer at lower frequencies since 
some of the radiation that would be scattered by an in- 
finite screen is now lost, causing a to be lower than the 
standard Kolmogorov value. Indeed, towards these direc- 
tions at low Galactic latitudes one encounters numerous 
HII regions embedded in the dense thin disk, so that the 
probability to find such multiple screens along the LOS is 
quite high. 

In order to get a better understanding of the scatter- 
ing properties in the ISM we need measurements of a for 
more LOSs. In particular the scattering properties of in- 
termediate DM pulsars (100 < DM < 400 pc cm" 3 ) 
have not been studied yet. In an effort to find a for in- 
termediate DM pulsars we have currently launched an 
observational program using the Gaint Metrewave Radio 
Telescope (GMRT) in Pune, India. In this paper we report 
results for the first phase of these observations. 

2. Observations and Data analysis 

The following selection criteria were used to find an ade- 
quate sample of intermediate DM pulsars for our GMRT 



observations. Firstly pulsars with flux density > 10 mJy 
at 400 MHz were chosen from the Taylor et al. (1993; 
updated version 1995) catalogue. We also ensured that 
the expected width of the pulse, i.e. the intrinsic width 
along with the pulse broadened width, at the lowest fre- 
quency (243 MHz) is smaller than 80% of the pulse pe- 
riod. Using an intrinsic pulse width of 5% of the period at 
~5 GHz we estimated the width at lower frequencies ac- 
cording to a power-law index of A 025 (e.g. Thorsett 1991), 
where A is the wavelength in meters. We found the ex- 
pected pulse broadening at each frequency by applying 
the empirical relation T d (ms) = 4.2 x 10~ 5 DM 16 x (1 + 
3.1 x 10" 5 DM 3 ) A 44 (Ramachandran et al. 1997). The 
final estimated width was found by adding the two widths 
in quadrature. Applying the above selection criteria we 
ended up with a sample of 33 pulsars. Here we report 
results of the first phase of GMRT observations of nine 
pulsars. 

The observations were carried out with the GMRT in 
February 2002 using three bands around center frequen- 
cies at 243, 325 and 610 MHz. The GMRT has a 'Y' shaped 
hybrid configuration of antennas with 14 antennas placed 
randomly in a compact central array of 1 km by 1 km, 
and the remaining 16 antennas distributed along the three 
arms of the 'Y' (Swarup et al. 1997). The GMRT was in 
its commissioning phase during our observations, and due 
to various maintenance activities not all of the 30 anten- 
nas were available for observations. The observations were 
carried out with typically 20 to 25 antennas. Dual circu- 
lar polarization signals from all the selected antennas were 
incoherently added (i.e., signals from each antenna were 
first detected and then added) in a 256 channel filter bank 
with a total bandwidth of 16 MHz. The summed signals 
were integrated to a time constant of 0.516 ms and were 
recorded for off-line analysis after adding the two polar- 
izations (see Gupta et al. 2000 for more details about the 
pulsar mode of operation of the GMRT). The filterbank 
characteristics result in dispersion smearing per channel 1 , 
£dm> for the pulsars observed of 5.7 ms < £dm < 14-5 ms 
at 243 MHz, of 2.4 ms < t DM < 6.1 ms at 325 MHz, and 
of 0.4 ms < t DM < 0.9 ms at 610 MHz. 

Typical observation times were between 10 and 45 min, 
depending on the observing frequency and flux density of 
the individual source. During off-line reduction the filter- 
bank signals were de-dispersed and summed. We analysed 
the signals for RFI and skipped those parts of the data 
with spikes larger than a factor of three of the rms of the 
data. The signals were folded with the topocentric pulse 
period to produce total power profiles. 

At 408 and 1408 MHz we used pulse profiles 2 observed 
with the 76 m Lovell telescope at Jodrell Bank, UK. Using 



1 £dm = 8.3 DM Aubw/i^ 3 fJ-s, where DM is the dispersion 
measure (pc cm" 3 ), v is the observing frequency (GHz) and 
Avbw the filterbank channel bandwidth (MHz). 

2 These profiles are publicly available from the European 
Pulsar Network (EPN) archive maintained by the MPIfR, 
Bonn, http://www.mpifr-bonn.mpg.de/div/pulsar/data/ 
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Lovell, 1408 MHz 



GMRT, 610 MHz 



Lovell, 408 MHz 



GMRT, 325 MHz 



GMRT, 243 MHz 




Longitude (deg) 

Fig. 1. Integrated pulse profiles and best-fit model pro- 
files for PSR B1831— 03 at different frequencies. The pro- 
files at 243, 325 and 610 MHz were observed with the 
GMRT, whereas the 408 and 1408 MHz profiles were taken 
from the EPN database (Lovell observations). The align- 
ment of the profiles for different frequencies was done with 
respect to the peak of the main pulse. 

both circular polarizations the signals were incoherently 
de-dispersed and added in filterbanks with varying band- 
widths. At 408 MHz 32 channels of 0.125 MHz bandwidth 
(2.4 ms < tnu < 6.0 ms) were used, and at 1408 MHz 32 
channels of 1 MHz each (0.5 ms < £dm < 1.2 ms) were 
employed. Details of the system can be found in Gould & 
Lyne (1998). 

In order to measure the pulse broadening time t<\ we 
used the method described in Paper I. Firstly, we con- 
struct a pulse template as a sum of Gaussian components 
fitted to an observed high frequency profile, where the 
scattering for the given pulsars is negligible. For our sam- 
ple, the 1408 MHz profiles were adequate. For each lower 
frequency, we then find r^ from the best fit of the model 
profile, which is the convolution of the template with the 
dispersion smearing and the adopted PBFs, to the ob- 
served profile. The exact functional form for the PBF of 
the ISM is not known. We thus analyse the fits for three 
trial PBFs, (1) the PBF for a thin screen (PBFi) and 
(2) for a uniformly distributed medium (PBF2) in an ISM 
with Gaussian density fluctuations, given by (Williamson 
1972, 1973): 

PBFr(t) = exp(-i/r d ) U(t) (3) 

PBF 2 (t) = (^r d 3 /8i 5 ) 1 / 2 exp(-^ 2 r d /4t)[/(t) , (4) 

where U(t) is the unit step function, i.e. U(t < 0) = 
0,U(t > 0) = 1. (3) The third PBF is characterized by 
density fluctuations with a Levy probability distribution 



function that has a power-law decay (Boldyrev & Gwinn 
2003) and an asymptotic form PBF 3 (t) = (i/r d ) _4/3 U(t) 
for i/rd 3> 1. As for the high DM pulsars presented in 
Paper I, we again find PBFi to be most appropriate to 
describe the observed scattering. In particular, the fits us- 
ing PBF2 and PBF3 cannot reproduce the long "scattering 
tails" observed at lower frequencies, resulting in \ 2 values 
that are larger by factors of 2 and more. Using the thin 
screen approximation we obtain best-fit values and uncer- 
tainties for Td from the \ 2 contours in the plane of Td and 
offset in phase. 

In Fig. Q] observed and best-fit model profiles for the 
observed frequencies are shown for PSR B1831— 03. The 
template is constructed from the 1408 MHz Lovell profile 
using two Gaussians. Note the high S/N ratio and qual- 
ity of the GMRT profiles proving that this new telescope 
is highly capable of pulsar observations at low radio fre- 
quencies. The best-fit model profiles describe the shape of 
the observed profiles in an excellent manner. At 610 MHz 
the small peak at the leading part of the profile was not 
observed, which, however, does not affect the Td measure- 
ment (see next paragraph). At 243 MHz the dispersion 
smearing at the leading part of the profile seems to be 
not adequately described by the model profile, resulting 
in a much steeper rise of the peak. We repeated the fit 
using artificially increased dispersion smearing functions 
and found that the effect on Td is well below its la uncer- 
tainty and therefore negligible. 

As noted in Paper I, intrinsic profile variations with 
frequency (see the 610 MHz profile of Fig. Q could in 
principle give rise to inaccurate estimation of pulse broad- 
ening times. A careful analysis of these effects on the mea- 
sured Td can be done using simulated pulse profiles with 
frequency evolution that are made subject to pulse broad- 
ening. As shown, resulting deviations of the measured Td 
values from the true ones are in fact very small and can be 
accounted for using increased error bars. Thus, we again 
quote conservative 3tr error bars for all scatter broadening 
times. 

Recently, another method to analyse pulse broadening 
related to the CLEAN alogrithm was proposed by Bhat 
et al. (2003). In their approach, the authors try to derive 
the intrinsic pulse shape at the observed frequency with- 
out using any knowledge of the pulse profile at another, 
higher frequency. They point out that utilizing a high fre- 
quency template can indeed lead to uncertainties due to 
the same unkown frequency evolution of the pulse profile 
that we try to simulate in our computations (see Paper I) . 
Whilst it is indeed more straightforward in their method 
to perform a deconvolution to recover the intrinsic pro- 
file, their alogrithm cannot always produce unique results, 
yielding strikingly different values and hence uncertain- 
ties, sometimes. This is demonstrated for PSR B1849+00 
which was also studied in Paper I. Applying PBFi and 
PBF2 (see Eq. [3] and QJ the authors obtain equally good 
fits for Td = 225 ± 14 ms and Td = 121 ± 6 ms, where a 
choice can only be made by making an assumption about 
the more likely intrinsic profile. A comparison of these val- 
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Table 1. Pulse Broadening Measurements 



PSR 


I 


b 


DM 


D 


T243 


T~325 


1"408 


7610 


Q 


P 


logC n 2 o 




(deg) 


(deg) 


(pc cm -3 ) 


(kpc) 


(ms) 


(ms) 


(ms) 


(ms) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


B1821-19 


12.3 


-3.1 


224 


4.7 




58 (23) 


22 (4) 


4.6 (0.2) 


40±° ;? 


4.0! -.4 


-1.1 


B1826-17 


14.6 


-3.4 


218 


4.7 




137 (37) 


41(6) 


7.5 (1.3) 


4.61^ 


o fi +0.4 
°' D -0 .2 


-0.9 


B1831-03 


27.7 


2.3 


236 


5.2 


208 (94) 


59 (4) 


21.3 (1.1) 


3.2 (0.5) 


4.5±° 1 


o fj+0.5 
O.D_ 2 


-1.2 


B1845-01 


31.3 


0.0 


159 


4.0 




262 (181) 


54 (8) 


17.0 (1.4) 


4.3±?:£ 


O Q + 1.9 

J.O_ 4 


-0.5 


B1859+03 


37.2 


-0.6 


401 


7.3 




241 (51) 


89 (9) 


12.6 (0.9) 


4 7 +o.3 


q c+0.2 


-1.0 


B1900+01 


35.7 


-2.0 


246 


3.4 


176 (58) 


34 (2) 


14.7 (1.3) 


< 1.0 


4.8±° -j 


3.4±°oJ 


-1.0 


B1920+21 


55.3 


2.9 


217 


7.6 


4.5 (1.5) 


<0.8 


< 1.7 


< 0.6 






-2.9 


B1933+16 


52.4 


-2.1 


159 


5.6 


4.6 (0.2) 


1.8 (0.1) 


< 1.9 


< 0.4 


3-4±° j 




-2.7 


B2002+31 


69.0 


0.0 


235 


7.5 


27 (13) 


5.3 (1.2) 


2.0 (1.4) 


< 0.7 


S.Oli'I 


q q+13 
°-°-0.S 


-2.3 



Note.- Cols. (2) and (3) give the Galactic longitude and latitude of each pulsar; col. (4) its Dispersion Measure (DM) and 
col. (5) its DM distance (from the Cordes & Lazio (2002) model, hereafter NE2001); cols. (6)-(9) the pulse broadening times 
with 3<r errors in parenthesis at 243, 325, 408 and 610 MHz. Upper limits are quoted with "<". Col. (10) gives the spectral 
index of pulse broadening and col. (11) the spectral index of density irregularities (using /3 = 2a /(a — 2)) for each pulsar with 
their la errors. Col. (12) gives the logarithm of the fluctuation strength Cn e . 



ues with our measurement of t& = 223 ± 24 ms as derived 
in Paper I shows that both methods result in consistent 
pulse broadening times for the case of the thin screen ap- 
proximation. This supports our findings that an exponen- 
tial decay is the most appropriate form to describe pulse 
broadening for intermediate and high DM pulsars. The 
example of PSR B1849+00 shows that extra, a priori in- 
formation (typically an idea of the expected pulse shape) 
is usually needed to obtain correct solutions for more com- 
plicated profiles which holds true for both the CLEAN al- 
gorithm as well as our approach. Given the apparent im- 
perfections of both methods, all derived values should be 
treated with considerable care, e.g. by reflecting the pos- 
sible systematic errors by increasing the error estimates 
correspondingly, as done in our study. It is comforting to 
note that for PSR B1849+00 the frequency dependence of 
Td, derived by Bhat et al. (2003, a — 3.5 ± 0.7), and us 
(Paper I, a = 2.8^o q) are consistent. Recent OH observa- 
tions toward PSR B1849+00 revealed absorption features 
that most likely originate from a small and dense molec- 
ular clump (Stanimirovic et al. 2003). Thus, the LOS to 
the pulsar probes complex material so that our findings of 
non-Kolmogorov frequency dependence of pulse broaden- 
ing is not surprising. 



3. Results and Discussion 

Table ^ summarises the multi- frequency measurements of 
pulse broadening times for our sample of pulsars. PSRs 
B1821-19 and B1859+03 were not observed at 243 MHz 
due to the limited telescope time. For PSRs B1826— 17 
and B1845— 01 we did not detect a pulse profile at this 
frequency, probably because the pulse is smeared due to 
high scattering. 



Fig.|Ushows the measured frequency dependence of Td. 
For all pulsars with more than one Td measurement, i.e. for 
all but PSR B1920+21, we calculated the spectral index of 
pulse broadening, a, from Monte-Carlo simulations as de- 
scribed in Paper I. For PSR B1933+16, we used values of 
r d = 67(13) ms at 110 MHz and r d = 25(4) ms at 160 MHz 
from the literature (Rickett 1977; Slee et al. 1980). The 
median a and its la errors are listed in Tabled Also listed 
are the spectral indices of the electron density spectrum, 
/3, which we calculated from the relation (3 — 2a/ (a — 2). 
Note that this relation is only valid for f3 < 4, i.e. for all 
our pulsars except PSR B1933+16. The average spectral 
index for seven pulsars (excluding PSRs B1920+21 and 
B1933+16) is (a) = 4.57 ± 0.09. Fig. El shows the spec- 
tral index of pulse broadening as a function of DM for our 
measurements as well as published data. All of our new 
spectral indices, except the one for PSR B1933+16, are 
consistent within their errors with the Kolmogorov value 
of 4.4, as found for most of the low DM pulsars as well. 
However, along several LOSs, as towards the Crab and 
Vela pulsar, PSR B1933+16 and the high DM pulsars, 
the Kolmogorov theory fails leading to flattening of the 
spectra of scatter broadening. The results for the Crab 
and Vela pulsar are not surprising, as the complex struc- 
ture of the surrounding supernova remnants can introduce 
a number of effects affecting the observed scattering prop- 
erties (e.g. Backer et al. 2000; Lyne et al. 2001). 

The spectral index of pulse broadening for the LOS 
to PSR B1933+16 is a = 3.4(2), i.e. significantly lower 
than the Kolmogorov value of 4.4. This measured flatten- 
ing of the spectrum has been observed before for high DM 
pulsars and can be explained by anomalous scattering at 
multiple scattering screens with finite extensions and/or 
varying scattering strength (see Paper I) . Our pulse broad- 
ening measurements for PSR B1933+16 along with pub- 
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B1821-19 B1B26-17 B1B31-03 B1845 T 01 B1859+03 



100 200 500 100 200 500 100 200 500 100 200 

Frequency (MHz) 



500 100 200 



Fig. 2. Pulse broadening times, Td, with their 3a errors 
as a function of observing frequency, v, for nine pulsars. 
Arrows denote measured upper limits on Td. The lines 
correspond to the linear fit of the form y — ~ax + K 
where y = log(rd) in ms and x = log(^) in MHz. The 
dashed lines are examples of the expected dependence due 
to a Kolmogorov spectrum, i.e. a — 4.4. The bottom right 
most panel shows a linear function with the derived aver- 
age spectral index (a) = 4.57. See text for further details. 



lished Td's at 110 and 160 MHz fit excellently to a power- 
law model with log(Td(msec)) = —3.4 log(i/)+8.9, where v 
is in MHz (as seen in Fig.EJ. It should be noted that com- 
bining measurements at different epochs can be affected by 
refractive scintillations which may alter Td between epochs 
and thereby alter values of a. However, our Td measure- 
ments being quasi simultaneous (separated by days) are 
unlikely to be affected by refractive scintillations and the 
spectrum derived by these values is in good agreement 
with the low frequency values obtained by Cordes et al. 
(1985) at an earlier epoch. The authors quote A^d values 
for this pulsar near closeby frequencies of 1.41, 1.42, and 
1.67 GHz as 0.125, 0.100, and 0.110 MHz, respectively. We 
extrapolated our Td spectrum to these frequencies, solved 
Eq. ^for C\ and found it to be 9.1, 7.1, and 4.5, respec- 
tively, i.e. much larger than unity. The A^d measurements 
could have been biased by refractive ISS (e.g. Gupta et 
al. 1994) leading to an underestimation of the true band- 
widths. The resulting smaller values of Td at these fre- 
quencies, however, would lead to a steeper spectrum or to 
higher C\ values making the discrepancy even worse 3 . Our 
measured Ci's are much greater than unity and therefore 
in contradiction to what is predicted by standard theories 
of the ISM. Lambert & Rickett (1999), for instance, tabu- 
lated possible values for different geometries and spectral 
models and found values ranging from 0.56 to 1.53, i.e. still 
much smaller than those for PSR B1933+16. However, if 



3 Note that there is a possibility that the nature of the PBF 
may evolve from the well established exponential form, in the 
regime where Az^d is the easier quantity to measure. However, 
there is no evidence of this from the main theoretical models 
(e.g. Lambert & Rickett 1999 and references therein). 
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Dispersion Measure (pc cm ) 

Fig. 3. Spectral index of pulse broadening, a, as a func- 
tion of dispersion measure, DM, for the current sample as 
well as for earlier measurements. The dashed line a — 4.4 
indicates the spectral index for a Kolmogorov spectrum. 
The dotted line a = 4.0 represents the spectral index for 
a Gaussian distribution of irregularities. 



we were to apply a Kolmogorov spectrum (a = —4.4) - 
even though it is a very poor fit to the data - the extrap- 
olated Td's for PSR B1933+16 would be much smaller at 
the corresponding frequencies, resulting in Ci's that are 
much smaller, eventually even smaller than unity. Thus, 
the LOS to PSR B1933+16 remains to be an unresolved 
puzzle. Careful observations of decorrelation bandwidths 
of this pulsar in the frequency range of 300—1400 MHz are 
needed to help understanding the properties of the ISM 
along this LOS. 

One of the conventional methods of examining the 
scattering properties along various LOSs is to investigate 
the relation between the measured pulse broadening and 
the DM. Considering DM as proportional to the distance 
D to the pulsar, the power law model of density irregu- 
larities predicts r d oc (C£j 2 /0»-2) DM' 3 /^- 2 ' l/ -2/9/(/S-2) j 

which reduces to r d oc (C 2 ^) 1 - 2 DM 22 v~ AA for the Kolmo- 
gorov spectrum with j3 — 11/3 (Romani et al. 1986). 
Sutton (1971) first noted that, for DM < 20, Td increases 
roughly as DM 2 , but for 20 <, DM < 400 the rela- 
tion steepens approaching DM . With many more avail- 
able measurements of Td (and A^d) this relation has been 
revisited (e.g. Rickett 1977; Ramachandran et al. 1997). 
However, the basic feature of flatter and steeper slopes for 
low and high DM pulsars seems to hold well. In most of 
these studies the measured Td or Afd are scaled to a refer- 
ence frequency assuming a a priori frequency dependence 
of v iA . In Fig. 01 we show a plot of Td versus DM at 400 
MHz for all the pulsars where the frequency scaling is now 
measured, i.e. all pulsars from Fig. [31 We have fitted the 
empirical function Td (ms) = ADM 7 (1 + BDM^), as was 
suggested by Ramachandran et al. (1997). Here the term 
(1 + BDM^) should provide a useful description for the 
apparent mean dependence of the level of turbulence on 
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DM. In order to assess the deviation from the Kolmogorov 
value we fixed 7 = 2.2, and derived a best fit of the form 



r d (ms) = 2.5 x 1CT 6 DM 22 (1 + 1.34 x 1CT 4 DM 



2.3\ 



(5) 



which is shown by solid line in Fig. 0] Comparing our re- 
sults with the fit obtained by Ramachandran et al. (1997), 
their best-fit curve is slightly shifted to higher r^ values, 
most probably caused by the large scatter of data points. 
The exponent £ is somewhat smaller than their value of 
2.5, leading to a flattening of the curve at the highest 
DMs. As the high DM pulsars have frequency scaling laws 
which are significantly smaller than 4.4 (see Paper I), the 
extrapolation of their spectrum down to 400 MHz leads 
to lower Td values than for the case of a Kolmogorov spec- 
trum. This explains the flattening of the Td-DM curve at 
the high DM end. However, the large scatter in the Td 
values and small number statistics do not allow a tighter 
constraint of £. The different slopes for low and high DM 
pulsars have been studied in detail by several investiga- 
tors (e.g. Rickett 1977; Cordes et al. 1985; Rickett 1990). 
An extra bias could be introduced by a wrong conversion 
of Az^d to Td for the low DM pulsars. Indeed, if C\ ~ 1 
(see Eq. ^| does not hold, a comparison of the /S.v& and 
Td measurements is not possible. However, if C\ differs 
from unity by a constant amount, all points are shifted in 
the Td versus DM plot to one direction by a fixed factor. 
Only a systematic change in C\ with DM could change 
the slope of the curve itself. However, as was discussed by 
Rickett (1977), theoretically there is very little flexibility 
in changing C\ from values close to unity. The only ex- 
perimental confirmation for this became available for the 
Vela pulsar by Backer (1974), who found C\ = 1.07 by ex- 
trapolating several Td and Ap^ measurements at a given 
frequency. For PSR B1933+16, however, as discussed ear- 
lier, it appears that there is some evidence for C\ values 
greater than unity. 

Alternatively, if we accept C\ ~ 1, one can in prin- 
ciple explain the steeping in the Td versus DM relation 
by proposing a breakdown of the homogeneity condition 
in the medium (Rickett 1970). A usual way to check this 
condition is to analyse the scattering strength C 2 along 
the LOS. For a homogeneous medium with a Kolmogorov 
spectrum, C 2 can be found by using the definition from 
Cordes (1986), as 



CL (m- 



-20/3\ 



0.002 i/iVSp-ii/e ^jB/a 



(6) 



with v in GHz, D in kpc and Ai/d in MHz. We used Eq.El 
to calculate C 2 for the pulsars from this paper and Paper 
I, by averaging over all frequencies. Our Td values were 
converted to A^d using C\ — 1.16, the value appropriate 
for a uniform medium and having a Kolmogorov spec- 
trum (Cordes & Lazio 2002). The distances were obtained 
from the Cordes & Lazio (2002) model of the Galactic 
electron density distribution, hereafter NE2001. C 2 for 
other pulsars were used from Johnston et al. (1998) and 
Cordes et al. (1985), but scaled appropriately for the new 
NE2001 model distances. In Fig. we plot the scattering- 




Fig. 4. Pulse broadening times, Td, versus dispersion mea- 
sure, DM, for 27 pulsars as plotted in Figure Filled 
circles corresponds to decorrelation bandwidth measure- 
ments and open circles are pulse broadening measure- 
ments. The solid line corresponds to the best fit given 
by Eq. [SJ See text for further details. 



strength C 2 as a function of distance for all the 27 pul- 
sars with known frequency scaling law. As clearly seen in 
the top panel, C 2 increases with distance, a result which 
is consistent with that of Cordes et al. (1985). For a ho- 
mogeneous medium one expects to find log(C 2 ) near its 
canonical value of —3.5 (e.g. Johnston et al. 1998), which 
is only true for the nearby pulsars. Instead, for most of 
the sources with D ^ 3 kpc, we find values much larger 
than —3.5 indicating that we sample highly dense regions 
with enhanced scattering along these LOSs. As seen in 
the bottom panel of Fig. [21 the LOSs with enhanced scat- 
tering strength are concentrated in the inner regions of 
the Galaxy where the majority of HII regions can be 
found. These LOSs additionally show flatter spectra of 
pulse broadening indicating multiple scattering screens of 
anisotropic irregularities (Paper I). 



4. Summary & Conclusions 

We have presented measurements of pulse broadening 
times for a sample of intermediate DM pulsars. The ob- 
tained results bridge the gap in the available literature 
from low DM to high DM pulsars. For both the high DM 
pulsars presented in Paper I and the current sample, we 
find that the PBF for a gaussian distribution of irregular- 
ities and applying the thin screen approximation is most 
appropriate to describe the observed scattering. While 
there are significant deviations from the expected fre- 
quency scaling for high DM pulsars (Paper I), our current 
sample is, apart from the exception for PSR B1933+16, 
consistent with a Kolmogorov spectrum. Therefore, we 
conclude that LOSs to pulsars with Si DM is 300 pc 
cm" 3 show pulse broadening that is consistent with a 
Kolmogorov spectrum of electron density irregularities 
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Fig. 5. The top panel shows scattering strength log(C^ e ) 
versus distance for 27 pulsars as in Figure [3J The dashed 
line corresponds to —3.5, the canonical value for a homoge- 
neous Kolmogorov medium. The bottom panel shows the 
LOSs for these pulsars in the Galaxy, where the size of the 
circles corresponds to the magnitude of C^ e . The filled cir- 
cles indicate decorrelation bandwidth measurements and 
the open circles are pulse broadening measurements. See 
text for further details. 



(except for a few pulsars with complicated LOSs, e.g. 
the Crab and Vela pulsars, PSR B1933+16). At around 
DM= 300 pc cm -3 a change in the spectral index of pulse 
broadening is observed leading to a flattening of the spec- 
tra. We think that this change is related to a change of the 
Galactic material in the inner region of the Galaxy. Future 
multi- frequency observations are highly desirable to probe 
the DM range of 250 — 400 pc cm~ 3 and to determine the 
exact transition point. 

Standard theory appears still to be challenged by our 
result for the relationship connecting the pulse broadening 
time with the decorrelation bandwidth. Mostly only one of 
these quantities is measurable, while the other is computed 
using an assumed standard relation. This can obviously 
lead to systematic errors which need to be considered 
when combining corresponding data sets. Nevertheless, 
such procedure may be still unavoidable in order to in- 
crease the size of the studied sample, since the results ob- 
tained here and in Paper I underline the impression that 
the overall state of the ISM can only be determined in 
a statistical sense. Clearly, individual results are affected 
by the properties of certain LOSs. For this reason, special 
care is needed when interpreting the results obtained only 
for a small sample of pulsars. 
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